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Abstract 

A ( 7 -ary maximum distance separable (MDS) code C with length 
n, dimension k over an alphabet A of size q is a, set of codewords 
that are elements of A^, such that the Hamming distance between two 
distinct codewords in C is at least n — k + 1. Sets of mutually orthog¬ 
onal Latin squares of orders q <9, corresponding to two-dimensional 
g-ary MDS codes, and g-ary one-error-correcting MDS codes for g < 8 
have been classified in earlier studies. These results are used here 
to complete the classification of all 7-ary and 8 -ary MDS codes with 
d > 3 using a computer search. 


1 Introduction 

A q-ary code C of length n, and size M is a set of M elements, called code¬ 
words, of A", where A is an alphabet of size q. The minimum distance d of a 
code is the smallest Hamming distance between any two distinct codewords. 
A code with these parameters is called an (n, M, d)q code. If A is a finite 
held and U is a vector subspace, then C is called linear. A code that is not 
linear is called nonlinear. Codes that can be either linear or nonlinear are 
called unrestricted. 

In the unrestricted case, two codes are called equivalent if one can be 
obtained from the other by a permutation of coordinates followed by permu¬ 
tations of symbols at each coordinate separately. These operations preserve 
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the Hamming distances between codewords. We say that an (n, M, d)g code 
is unique if all codes with the same parameters are equivalent. 

An upper bound for the size of an (n, M, d)q code is the Singleton bound, 

M < 

Codes meeting this bound are called maximum distance separable (MDS), 
and k = n—d+1 is called the dimension of the MDS code. A q-aiy MDS code 
with length n and dimension k is called an (n, k)q MDS code. MDS codes 
have the property that, given any k coordinates, each fc-tuple of symbols 
from A occurs in the given coordinates in exactly one codeword. It has been 
conjectured (the MDS conjecture) that an MDS code with parameters n, k, 
and prime power q with 1 < A: < n — 1 exists if and only if n < q + 1, with 
the exception that when g is a power of two, MDS codes with n = q + 2, 
and k = 3oTk = q — 1 exist [21]. The conjecture has been proved for linear 
codes when g is a prime by Ball [Ij and when g is a power of a prime p and 
A: < 2p — 2 by Ball and De Beule [S] . 

MDS codes with d = 1 or A; = 1 are unique and they are called trivial; 
the first contains the whole space and the latter is a repetition code. The 
case k = 2 corresponds to sets of mutually orthogonal Latin squares, which 
have been classified for O' < 9 [7|. MDS codes with d = 2, dimension k and 
alphabet size q correspond to A:-dimensional Latin hypercubes of order q. For 
q = 2,3, they are trivially unique. Potapov and Krotov 123 give a recursive 
formula for the number of (A: + 1, k)^ MDS codes but do not classify them 
up to equivalence. With small k and q, {k + l,A:)q MDS codes have been 
classihed by McKay and Wanless [Hj. 

For q = 2, nontrivial MDS codes with d > 2 do not exist. For q = 3, 
the only nontrivial MDS code with d > 2 is the unique (4, 2)3 MDS code. 
Alderson [1] showed that the ( 6 , 8)4 and the (5, 8)4 MDS codes are unique. 
The nonexistence of pairs of mutually orthogonal Latin squares of order 6 
implies the nonexistence of nontrivial 6 -ary MDS codes with d > 8. For 
q = 5,7 all MDS codes with d > 8 , except the ( 4 , 2)7 codes, are equivalent 
to linear codes [lO]. For q = 5, they are unique, which follows from the 
uniqueness in terms of the notion of equivalence of linear codes | 6 |. For 
q = 7, S, the MDS codes with d = 8 were classified in [TOl |TT]. For q = 7, 
some classihcation results exist for codes with d > 8 in terms of different 
notions of equivalence for linear codes; see for example [9]. 

In this work, we classify all 7-ary 8 -ary MDS codes with n > 3, d > 3 
by a computer search to finish the classification of all 7-ary and 8 -ary MDS 
codes with d > 8 . In Section [2l we discuss basic properties of MDS codes 
and computational tools used in this work. Section [3] explains the computer 


2 


search used, and the results are given in Section HI Finally, in Section 0 we 
discuss the case d = 2 corresponding to Latin hypercubes for small q. 

2 Preliminaries 

2.1 MDS codes 

The operations maintaining equivalence of unrestricted codes form a group 
of order (g!)”n!, denoted by that acts on A^. The stabilizer of a code C 
under this action is called the automorphism group of the code. Ant((7). 

Let C be an (n, k)q MDS code. For a symbol i E A, we denote by Ci 
the {n — 1, k — l)i} MDS code obtained by removing the last coordinate and 
retaining the codewords that have symbol i in that coordinate, 

Ci = {(Ci, C2, . . . , Cn-l) : cEC,Cn = i}. 

Removing a symbol at a given position of each codeword of C yields 
an (n — l,k)q MDS code C. This operation is called puncturing. We say 
that C is an extension of C. An (n, k)q MDS code C is extendable if there 
exists an extension of C that is an {n + 1, k)q MDS code. There is a one-to- 
one correspondence between extensions of C (up to the position of the new 
coordinate) and labeled partitions of C into (n, k — l)q MDS codes: Uieyi 
where CH denotes the code obtained by adding the symbol i at the end of 
each codeword in (7*, is an extension of C if and only if the sets (7* form a 
partition of C into (n, k — 1 )^ MDS codes. 

It is known that if a linear code is extendable, it has an extension that is 
linear [2]. For hxed q and d and sufficiently large n, an extension of a linear 
MDS code is necessarily equivalent to a linear MDS code [3] . However, very 
little is known about extending nonlinear codes. 

2.2 Tools 

We reduce the problem of detecting code equivalence to the graph isomor¬ 
phism problem and use the software nauty [TB] for solving the instances. For 
an (n, M, d)q code (7, we dehne a colored graph as follows: The graph con¬ 
tains n copies of Kq, the complete graph of order q, named Fi,F 2 ,... ,F„, 
colored with color 1. For each i, the vertices in F, correspond to the elements 
in A. For each codeword c E C, there is an additional vertex, colored with 
color 2, which is adjacent to the vertex corresponding to the symbol c* in 
Fj for each i E {1,2,... ,n}. A graph isomorphism preserving the coloring 
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permutes the copies of Kq and the vertices in each of them separately, cor¬ 
responding to a permutation of coordinates and permutations of symbols in 
each coordinate separately in the code, respectively. Two codes are equiva¬ 
lent if and only if their corresponding graphs are isomorphic. Further, the 
automorphism group of a code corresponds to the automorphism group of 
the graph. The software nauty can be used to hnd the automorphism group 
of a graph and all isomorphisms between two graphs. For large graphs, we 
use nauty in the sparse mode with the random Schreier algorithm enabled. 

For a hnite set X and a family S of subsets of X, an exact cover is a 
subset S' of S that partitions X. We use the library libexact [8] for hnding 
all exact covers, given X and S. 

Finally, we use cliquer Cl to hnd all cliques of given size in a graph. 


3 Generation 

Given a set C of equivalence class representatives of (n, k)q MDS codes, every 
(n -|- 1, k)q MDS code is equivalent to a code that is an extension of a code in 
C. This reduces the problem of generating equivalence class representatives 
of {n+ 1, k)q MDS codes to the problem of hnding all partitions of all codes 
in C into (n, k — l)g MDS codes. For k = 2, this is precisely the method of 
extending a set of MOLS with a new Latin square by partitioning the Latin 
squares in the set into common transversals, used for example in [H |T31 HU 
[T^ : for discussion of this and similar problems, see for example [T 8 ] . 

Consider a partition of an (n, k)q MDS code C into (n, k — 1)^ MDS codes 
(7* for z e M, 

c=\Jc\ 

ieA 

We have 

c;, = U F 

ieA 

and 

= u Fj- 

j&A 

Finding all possible partitions of C into (n, k — 1)^ MDS codes can now 
be reduced to hnding all possible partitions {Cj}i^_A of Cj for each j and 
combining those partitions in all possible ways. Instead of solving the reduced 
problem using recursion, our method for hnding the partitions of (n, k)q MDS 
codes uses the results from hnding the partitions of (u —1, k — l)g MDS codes. 
In this work, we generate all g-ary MDS codes with k > 3, d > 4 starting 
from the MDS codes with k = 2, their partitions, and the codes with d = 3. 
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3.1 Algorithm 

As an initial step, we find all partitions of the eqnivalence class representa¬ 
tives of (n, 2)q MDS codes into (n, l)q MDS codes in the following way. Let 
Cn be the set of representatives of (n, 2)^ MDS codes. For a given n, we loop 
over all C G Cn+i and over all n -|- 1 pnnctured codes C of C. Now, C is an 
extension of C so it indnces a partition of C into (n, l)g MDS codes. We 
hnd the eqnivalence class representative C' G for which C' = C and hnd 
all g G Gn for which gC = C'. Considering how the codewords are mapped 
by g, we get partitions of C' into (n, l)q MDS codes. All partitions of each 
element in C„ are obtained in this way, as every snch partition corresponds 
to an (n -|- 1, 2)q MDS code C that is eqnivalent to some C G C„+i. 

The algorithm for hnding all partitions of an (n, k)q MDS code C into 
(n, k — l)q MDS codes consists of three parts. 

1. For each j G A, hnd all {n — l,k~ 2)g MDS codes that occur as a part 
in a partition of Cj into {n — l,k ~ 2)^ MDS codes. Because Cj = C, 
where (7 is a representative whose partitions into {n — l,k ~ 2)^ MDS 
codes are known, this can be solved by hnding & g E Gn for which 
Gj = gC and applying g directly to each partition of C. Denote the 
set of (n — 1, A; — 2)q MDS codes found in this step by Sj. 

2. Find all subsets D of G that are (n, k — 1)^ MDS codes for which 
Dj G Sj for each j. This is done by hnding all g-cliques in the q'-partite 
graph where, for each j, each element of Sj corresponds to a vertex 
in the jth part, and between two vertices in diherent parts there is an 
edge if the minimum distance between them is at least n — k + 1. 

3. Find all partitions of G into (n, k — 1)^ MDS codes. This is done 
by solving an exact cover problem, where q sets are selected from the 
family created in the previous phase such that each codeword in G is 
covered exactly once. 

We run the algorithm for one representative G from each equivalence class 
of (n, k)q MDS codes, and hnally perform isomorph rejection for all k)q 

MDS codes corresponding to the obtained partitions. 

3.2 Consistency check 

We check the consistency of each phase of the generation by double counting. 
The total number of {n -|- f,k)q MDS codes is the sum of the sizes of the 
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equivalence classes, 


E l^n+l| 

where Cn+i is the obtained set of equivalence class representatives of (n+ 1 , k)q 
MDS codes. On the other hand, the same count can be obtained as follows. 
Let N{C) be the number of partitions of an [n, k)q MDS code C into {n, k—l)q 
MDS codes. The number of (n + 1, k)q MDS codes is 


5^JV(C),!= 5^ 

C ceCr, 


Gn\N{C) , 
|Aut(C')| 


where the hrst sum is taken over all {n,k)q MDS codes, Cn is a set containing 
exactly one representative from each equivalence class of (n, k)q MDS codes, 
and N{C) is the number of partitions of C obtained by the algorithm. 


4 Results 

The algorithm was run for q = 7 and q = 8, for each k = 3,4,..., g — 1, 
and for each n starting from n = k + 2 and increasing n until no more codes 
were found. The numbers of inequivalent extendable codes are shown in 
Tables [T] and [3l for g = 7, 8, respectively, and the numbers of equivalence 
classes are shown in Tables [2] and 01 In the tables, the numbers marked 
with * are obtained in this work using the method described above. The 
previously known numbers are obtained from the classihcation of one-error- 
correcting MDS codes and sets of MOLS, and in addition for g = 7 the MDS 
conjecture gives an upper bound for n, and the nonexistence of (10, Sjs MDS 
codes follows form the nonexistence of (9, 2)3 MDS codes. The computations 
lasted a few minutes for g = 7 and 21 hours for g = 8 ; 18 hours were needed 
for extending the (9, 7)3 MDS codes. 

By comparing the equivalence class representatives with known linear 
codes, we see that all 8 -ary MDS codes with d>5oTk> 4, (1 = 4 are 
equivalent to a linear code. 


5 MDS codes with d = 2 

The preceding part of this work completes the classification of all MDS codes 
with minimum distance at least 3 for alphabet sizes at most 8 . In this section, 
we discuss the remaining case of minimum distance 2. We use a simple 
construction to give a lower bound for the number of equivalence classes 
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n\k 

2 3 4 5 6 7 

3 

6 

4 

2 1 

5 

1 *1 1 

6 

1 *1 *1 1 

7 

1 *1 *1 *1 1 

00 

0 0 0 0 0 0 


Table 1: Number of inequivalent extendable nontrivial {n,k)i MDS codes 


n\k 

2 3 4 5 6 7 

3 

147 

4 

7 

5 

1 1 

6 

1 *3 1 

7 

1 *1 *1 1 

8 

1 *1 *1 *1 1 

9 

0 0 0 0 0 0 


Table 2: Number of equivalence classes of nontrivial {n,k)i MDS codes 


n\k 

2 

3 

4 

5 

6 

7 

3 

2 024 






4 

38 

4470 





5 

1 

*44 

36 




6 

1 

*2 

*1 

12 



7 

1 

*2 

*1 


7 


8 

1 

*2 

*1 

*1 

*1 

4 

9 

0 

*2 

*0 

*0 

*0 

*1 

10 


0 




*0 


Table 3: Number of inequivalent extendable nontrivial (n, k)s MDS codes 
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n\k 

2 3 4 5 6 7 8 

3 

283 657 

4 

2165 

5 

39 12 484 

6 

1 *39 14 

7 

1 *2 *2 8 

8 

1 *2 *1 *2 4 

9 

1 *2 *1 *1 *2 4 

10 

0 *1 *0 *0 *0 *1 0 

11 

0 *0 


Table 4: Number of equivalence classes of nontrivial (n, k)s MDS codes 


of (n, n — l)q MDS codes, or equivalently the paratopy classes of (n — 1)- 
dimensional Latin hypercubes of order q, for small q. 

We denote by W the number of (n, n — 1)^ MDS codes and by Mn the 
number of equivalence classes of {n,n — 1)^ MDS codes. Because each class 
contains at most |Gn| = nl^ql)^ elements, we get 


Mn > 


Nn 

n\{ql)^ 


( 1 ) 


For small values of n and q, we hnd lower bounds for and by the 
following construction. This construction resembles Construction X4 given 
in m Chapter 18.7]. We denote by CC the direct sum of codes C and C 
of lengths n and n\ respectively. 


CD = {(ci, C 2 ,..., c„, c'l, 4,..., 40 : c e C, c' e C"}. 

Theorem 1. Let C and C he (ni,ni — 1)^ and {n 2 ,n 2 — l)g MDS codes, 
respectively. Then 

C" = IJ CiC[ 

i&A 

is an (ni + ^2 — 2, ni + n 2 — 3)^ MDS code. 

Proof. The code C" contains qq'^i-‘^q'^ 2-2 _ gni+n2-3 codewords of length 
Ui + 712 — 2. Because the minimum distance of each Ci and C' is 2, the 
minimum distance of each CiC[ is at least 2. For distinct i,i, any codeword 
in CiC[ differs from any codeword in CjC'^ in at least one position in the hrst 
771 — 1 coordinates and in at least one position in the last 772 — 1 coordinates. 
Therefore, C has minimum distance 2. □ 
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n\q 

4 

5 

6 

7 

8 

3 

2 

2 

12 

147 

283 657 

4 

5 

15 

264 248 

> 4.8 X 10^ 

> 4.6 X lO^s 

5 

26 

86 

> 4.8 X 10^ 

> 2.3 X 10^3 

> 6.2 X 1025 

6 

4 785 

3102 

> 1.4 X 10^3 

> 9.5 X 10^® 

> 6.9 X 1035 

7 

> 3.5 X 10^ 

> 1.5 X 10^ 

> 3.1 X 10^^ 

> 3.3 X 1024 

> 6.6 X 1045 


Table 5: Number of equivalence classes of (n, n — l)g MDS codes 


Because the construction in Theorem [T] maps exactly q\ pairs of {C, 
to the same code, this yields a lower bound for the number of {n, n — 
MDS codes, 


Nn > max 

n' 


^n' ^n—n'+2 

q\ 


C) 

1 ). 

( 2 ) 


This lower bound along with the known numbers of (n, n — l)q MDS codes 
and dl]) allows us to hnd lower bounds for the number of equivalence classes 
of (n, n — l)g MDS codes for small n. 

Table [5] lists the known numbers of equivalence classes of (n, n — l)q MDS 
codes for 4 < g < 8 given in [16] and lower bounds obtained by ([I]) and o 
for small values of n. Since the lower bounds for the unknown values are 
rather large, except for a few entries, obtaining further classihcation results 
by explicitly constructing representatives from each equivalence class is not 
feasible in general. 

Finally, we note that 121 gives a double exponential lower bound for Nn 
when g > 4: 


logs Nn > 


nr. 

^ (g-3)(g-l) 


(n-l)/2 


for even g, 
for odd g. 


From dll) , it follows that that the number of equivalence classes for g > 4 
also has a lower bound that is double exponential in n. 
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